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When a carbon nanotube is truncated with a certain type of edges, boundary states localized near 
the edges appear at the Fermi level. Starting from lattice models, low-energy effective theories are 
constructed which describe electron correlation effects on the boundary states. We then focus on 
a thin metallic carbon nanotube which supports one or two boundary states and discuss physical 
consequences of the interaction between the boundary states and bulk collective excitations. By 
the renormalization group analyses together with the open boundary bosonization, we show that 
the repulsive bulk interactions suppress the charge fluctuations at boundaries and assist the spin 
polarization. 
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I. INTRODUCTION 

A single-wall carbon nanotube (CNT) is a fascinat- 
ing quasi-one-dimcnsional (quasi-lD) nano-scale mate- 
rial, which is a graphite sheet wrapped into a cylindrical 
form. Its electronic structure is basically well described 
by a one-electron tight-binding model with a single ir or- 
bital per atom. Depending on its geometrical shape, a 
large variety of electronic structures are realized. Es- 
pecially, it can be either a metal or a semiconductor 
depending on how the graphite sheet is wrappediS*^ 
The way of wrapping is specified by a chiral vector 
(N,M),N,M G Z. When TV - M = mod 3 is sat- 
isfied, a CNT is metallic, while it is gapped otherwise. 

An interesting consequence of its rich electronic band 
structure is the existence of the boundary states when 
the system possesses boundaries. For a CNT with zigzag 
or bearded edges, there appear states localized at the 
boundaries for specific values of the wave number along 
the boundaries^ It is a hallmark of the phase degree of 
freedom specific to quantum mechanical systems^ The 
existence of such boundary states raises an interesting 
question as to what kind of physical consequences they 
lead. For example, the electronic and magnetic prop- 
erties of nanographite in magnetic field^ or electronic 
transport through nanographite ribbon junctions 7 ^ were 
theoretically investigated. Furthermore, in the presence 
of electron-electron or electron-phonon interactions, the 
boundary states might trigger an instability as they form 
a flatband and a sharp peak in density of states at the 
Fermi energy for a 2D sheet geometry. Indeed, spin po- 
larization induced by the boundary states^*i2*ii*i£, or 
coupling with lattice distortions^*^ has been studied for 
a graphite sheet by several authors. The effects of ID 
low-lying excitations localized at the boundaries was also 
discussed for ribbon geometry^ 

The effects of bulk electron correlations have been 
extensively investigated for metallic CNT's without 
boundariesi 15 i 16 i 17 i 18 i 19 i 20 i 21 i 22 It is claimed that the most 
important forward scattering part of the Coulomb inter- 



action is well accounted for by the Tomonaga-Luttingcr 
(TL) liquid picture 2 ^, where low-lying excitations are not 
of Fermi liquid type, but bosonic collective excitations. 
Behaviors specific to TL liquid, such as a characteris- 
tic temperature dependence of conductance, have been 
indeed observed in recent experiments^ In TL liquid, 
boundary critical phenomena are known to be drastically 
different from the conventional Fermi liquid case when 
the system possesses boundariesi 2 ^* 2 ^**^* 2 ^* 2 ^ Then, we 
expect to see interesting phenomena for a thin metallic 
CNT with boundaries. The anomalous boundary physics 
in a metallic CNT within the TL liquid picture such as 
tunneling density of states^ 6 -, Freidel oscillation 2 ^, or lo- 
cal density of states^ - has been investigated previously, 
but without boundary states. 

The purpose of the present paper is to discuss elec- 
tron correlation effects for CNT's with edges that sup- 
ports boundary states. We consider (TV, —TV) CNT's with 
zigzag and bearded edges, for which boundary states ap- 
pear at the Fermi level for some values of the wave num- 
ber along the edges. Starting from lattice models with 
the Coulomb or the Hubbard interaction, we first estab- 
lish low-energy effective theories that describe correlation 
effects at boundaries. We then focus on a thin metallic 
CNT, where the boundary states interact with the collec- 
tive bulk excitations. By the renormalization group (RG) 
analyses together with the open boundary bosonization, 
we discuss the cases where one or two boundary states 
exist. 

For the case of two boundary states, we are especially 
interested in whether or not the boundary states exhibit 
spin polarization in the presence of gapless bulk excita- 
tions. The possibility of spin polarization was previously 
discussed for 2D geometry, i.e., in the limit of the in- 
finite tube radius TV — > +oo, by mean field theory^ or 
density-functional theory with local spin density approx- 
imation (LSDA)i°i. However, these treatments of elec- 
tron correlations can overestimate the magnetic instabil- 
ity. Also, when the system is wrapped into a ID cylinder, 
one needs further justification. A density matrix rcnor- 
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FIG. 1: A (TV, -TV) CNT with (a) zigzag and (b) bearded 
edges. Two-dimensional graphite sheets are wrapped around 
the x axis. The ovals indicate how to form a spinor c = 
(c. , Co), and dotted squares show a choice of unit cell for each 
types of edges, (c) Energy spectra for a (TV, —TV) CNT with 
zigzag edges for TV = 6, 9. (t = 1) Shaded regions represent 
bulk energy spectra for TV — > do. Allowed wave numbers 
along the boundaries for TV = 6, 9 are shown by vertical lines. 
Boundary states are denoted by x . 



malization (DMRG) study was also performedii which 
reports spin polarization at boundaries, but only for a 
semi-conducting CNT. When the bulk part of the sys- 
tem is gapless, it is not clear if the boundary states show 
spin polarization since the total spin (and charge) car- 
ried by the boundary states can dissipate into a bulk part 
of the system through the electron interactions between 
boundary states and bulk gapless collective modes. 



The paper is organized as follows: section[H]is devoted 
to a construction of low-energy effective theories which 
account for zero-energy boundary states. We recall some 
known facts on boundary states in section Til Al and the 
Coulomb or the Hubbard interaction is projected to the 
low-energy sector of the Hilbert space spanned by the 
boundary states and gapless bulk modes in section Til Bl 
In section III CI we perform open boundary bosoniza- 
tion to account for the bulk repulsive interactions and 
discuss its consequences on boundary physics. Finally, 
we present RG analyses for the case with one boundary 
states (section lTTT|) and with two boundary states (section 
II V|) . We conclude in section Ivl 



II. ZERO ENERGY BOUNDARY STATES AND 
INTERACTIONS 

A. Zero energy boundary states 



1. Zigzag edges 

Let us first consider a (TV, —TV) CNT with zigzag 
edges. Our starting point is the single particle tight- 
binding Hamiltonian defined on the honeycomb lattice 
[Fig. da)]: 

ff ki „ = ^2[c\i,j)T_c(i,j + l) + h.c. 



+J(i,j)T + c(i-l,j + l) + h.c. 
+J(i,j)T c(i,jj\, 



(2.1) 



where c = *(c.,c ), is a spinor made of electron an- 
nihilation operators cy defined on different sublatticcs 
• ,o. Coordinates of the spinors are labeled by the unit 
cell to which they belong and the location within the 
unit cell. Unit cells are chosen so as to be compatible 
with the shape of the edges, and are labeled by their y 
coordinate y = jao, j = 1, ••• , N y , where N y = TV is 
the total number of sites along the y axis and ao the 
lattice constant. Spinors are located at X — \/3aoi, 
i = 1, ■ ■ ■ , Nx within a unit cell with TVx being the total 
number of sites along the X axis. Hopping matrix ele- 
ments for the zigzag case are given by T + = (— t) 1 



T 



(-*) 





1 



and T = (-/) 



1 

1 







where t is 



the hopping integral. As specified by the chiral vec- 
tor, the periodic boundary condition is imposed along 
the edges, which renders the wave numbers along the 
edges quantized, k v aa = 2mit/N y ,m £ Z. The band 
structure is then composed of a set of ID modes, each 
of which is characterized by the wave number. Per- 
forming a Fourier transformation along the y axis as 



-- 



e yV Ck y (i), -Hki,, is decomposed as 



N x 

H(k y ) = ^2[cl y (i)V x c ky (i + l)+h, 
1=1 



c L(*)VoC fe (i) 



(2.2) 



- = T_e- ik y ao and V = T + T_e+' lk y a ° + 



where V- 

T+e -ifc yQ0 Thig 1D Hamiltonian can be seen as a 
ID chain with alternating hoppingsi^i As for bulk elec- 
tron states, the energy spectrum is gapless for two val- 
ues of k y , k y ao = ±fcoao = ±27r/3, whereas it is 
gapped for other k y . This is best seen by performing a 
gauge transformation c,(i) — > (— } l e +lkya °/ 4 c,(i), c (i) — > 
(— )' l e~ lkya °/ 4 c (i), which transforms the Hamiltonian as 
V x -> (- e -Mr fc »°°):r_ and V -> 2 cos(fcj / a /2)T . Then, 
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the alternation can be completely removed for k y = 
±fc , where the Hamiltonian is reduced to H(±ko) = 

-tYsiti [ c ±.i c ±,i+i + h - c ] b y defining c± ko ,.(i) =■ 
c±,i=2i-i, c±k ,o(i) =: c±j=2i, and 2N X =■ N x . 

When truncated with edges, the above one-parameter 
family of Hamiltonians supports zero-energy boundary 
states for some specific values of fc y r^ as a manifes- 
tation of a nontrivial bulk band structure. For zigzag 
edges, they appear for — ir < k y ao < —2ir/3,+2ir/3 < 
kyCio < +7r, within a bulk energy gap at the Fermi en- 
ergy. The localization length of the boundary states 
continuously increases as we vary k y ao from +7r(— 7r) to 
+27r/3(— 2tt/3). For k y ao = ±tt, the boundary states 
is completely localized at the boundaries i = 1 and Nx, 
while the localization length for k y a^ = ±27r/3 is infinite, 
at which the branch of the localized boundary states con- 
tinuously merges with the extended bulk states. If we ex- 
plicitly construct wave functions for the boundary states, 
they are given by 

/*,,.(*) <x (-2e i3fc ^/ 2 cosfc,a /2)^ 1 , f kyi0 (i) = 

(2.3) 

(— 7r < k y cio < — 2ir/3, +2ir/3 < k y ao < +ir), after the 
gauge transformation. 



2. Bearded edges 

We turn to a CNT with bearded edges. Atomic con- 
figurations realizing bearded edges for 7r electrons were 
recently proposed based on an ab initio calculation^! 
We choose a different unit cell from that for the zigzag 
case as indicated in Fig. Cfb). (X is now equal to 
2v / 3ao*- i = 1, ■ • • , Nx-) We obtain a set of Hamiltonian 
parametrized by k y as in the zigzag case with the hopping 
matrices Vx and Vq being given by V x = T_+T_e~ lkyaa , 
and Vq = Tq. Again, the bulk energy spectrum is gapless 
for k y = ±fco, for which the Hamiltonian can be trans- 
formed to a simple form iJ(±fco) = — t X)(=i [ c ± i c ± i+i + 
h.c] by a gauge transformation c.(i) — ► e +lkya ° l / 2 c,(i), 
c (i) — ► e +lkya0l / 2 c o (i). On the other hand, boundary 
states appear for different values of k y from the zigzag 
case, — 27r/3 < k y ao < +27r/3. The wave function of the 
boundary states can be explicitly constructed as 



fk y ,»(i) oc ( - 2cosfc y a /2) 



-0-1) 



f ky , (i) = (2.4) 



(— 27r/3 < k y ao < +27r/3), after the gauge transforma- 
tions. 



B. Interactions 

As the boundary states with different k y are all degen- 
erate at the Fermi energy, electron interactions have pro- 
nounced effects. In what follows, we will construct low- 
energy effective theories describing correlation effects on 



the boundary states. The interacting part of the Hamil- 
tonian is written as 

H-int = 2 VRR';pp';ss' 

RR' ;pp' :ss' 

xcl s (R)cl, s ,(R')c p , s ,(R')c ps (R), (2.5) 

where R = (i,j) runs over the lattice cites in the tight- 
binding model, p = •/o represents a sublattice, and s =f 
/ I is a spin index. As for Vrh' - pp ' - ss i , we take either the 
Hubbard interaction 

Vrr'-pp'- SS ' = USrriSpp'6s,-s' (2.6) 
or the unscreened Coulomb interaction 

e 2 /n 



VRB: 



' ;pp';ss' 



(2.7) 

where n is an effective dielectric constant of the system, 
x, = X, x = X + V3ao/2 (zigzag), \/3a Q (bearded), R 
is the tube radius, and r z ~ a characterizes the radius 
of p z orbital, serving as a short-distance cutoff. In Ref. 
ITU k is estimated to ~ 1.4, while r z is determined to be 
0.526 x a in Ref. 0, from the requirement that the on- 
site interaction in the original tight-binding model corre- 
sponds to the difference between the ionization potential 
and electron affinity of sp 2 hybridized carbon. Since the 
Coulomb interaction is unscreened in CNT's, the Hub- 
bard interaction is less realistic. However, as we will see, 
their differences are small for the correlation effects at 
boundaries, due to the special properties of the boundary 
states, although the long-range nature of the Coulomb 
interaction has fundamental effects for the bulk electron 
states. 

Focusing on a low-energy sector, we construct a effec- 
tive theory with boundary states at zero energy When 
a CNT is metallic, two ID gapless modes are also in- 
cluded, while we drop all the operators belonging to 
gapped bands: 



-{hi) 



y k 

l=^e i ^c Q ( J ) + -l=6( J , J ), 



'» a=± 



(2.8) 



where c±(i) := c±k (i), b is a linear combination of 
boundary states, b(i,j) = J2k y e ' ky V Sk y («) e ky , with 
f ky (i) being an eigen wave function localized at bound- 
aries, e\ and e k are creation and annihilation operators 

for a boundary state, and means the summation is 
restricted to the wave numbers for which a boundary 
state appears. We project the interaction H iat into the 
reduced Hilbcrt space spanned by the two gapless modes 
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c±, c ± and the boundary states e£ , e k . Substituting 
Eq. (|2.8|l . H iDt is decomposed into a number of terms, 
each of which is a four-fermion interaction made of ci- 
ther Cj_, c ± , e£ , or e k . However, not all of them give 
rise to a contribution due to the momentum conservation 
along the edges. For example, a four-fermion interaction 
composed of three bulk electrons and one edge electron 
(e.g., c\ rS c_ s ,c_ s ,b s ,) does not appear for the zigzag cases 
since the momentum carried by the bulk electron is ei- 
ther +2tt/3 or — 27r/3, while that carried by the edge elec- 
trons is never equal to ±2-7r/3, and hence a momentum 
mismatch occurs. It should be also noted that boundary 
states are nonvanishing only on one of sublattices p = • 
as seen from Eqs. <|2.3ll an d lj23}- Due to this special 
property of boundary states and also since the boundary 
states are localized at boundaries, the Coulomb interac- 
tion and the Hubbard interaction are not quite differ- 
ent. So we restricted ourselves to the Hubbard inter- 
action for a while as it is simpler. We write the pro- 
jected interaction as H iat — > + H° d t s °, where H^ lk is 
solely composed of bulk electron operators, while H^ t sc 
includes boundary states. We further decompose H* d f 
as H' A r = H ed r 4 + H° d r 3 + H° d r 2 + # od f BCl . Firs, H° d r 4 

mt nit int mt 1 nit " nit 

consists solely of boundary states, 



where 
J = 



/ , "-us 2 

P ^ t c as c as> 



s= E iajM 



-■kyS' 1 



Pe = E \fky\ 2e l v s e k 

ky ,S 



A = 



c +l c -1 +c-ic +T 



ky ky~\ ky J, 



(2.11) 



We omit sublattice indices and simply write c± = c± m , 
b = b, henceforth, since p = o does not appear. Finally, 

H° d t eo1 and H° d t E ° 3 are given by 



h: 



u 



E [ c It c -t { c ii b -i + b \ c +i) + h - 



h3 

V- 



c\.b ,+b\c + Ac\ l c_ l + h.c. 



v ■ 



1,3 



+ (cL|6| + 6|c + |)6j6| + h.c. 



(2.12) 



E S(-ki + k 2 -k 3 + k 4 ) 



y k!,k2,k 3 ,ki 



Xe k 1 ^ e k 



:T e Ll e fc 4 i E & (O/fca (*')/k. W/fc, (*). 

(2.9) 



where we simply wrote fk,p=» as //. and <5(fc) is a lat- 
tice delta function; i.e., S(k) is equal to unity only when 
k is integral multiple of 2ir, while it is vanishing other- 
wise. When a CNT is insulating, it is enough to consider 
only -ff'jft 8 " 4 for boundary physics. Then, the problem is 
reduced to how the degeneracy between the boundary 
states with quenched kinetic energy is lifted by iJ.^ 8 " 4 , 
as in the flatband magnetism or the fractional quantum 
Hall effect. Note also that the strength of the Hubbard 
interaction is reduced by the factor 1/N y , as the wave 
functions of the boundary states are extended over the 
circumference of the tube. 

The part composed of two edge operators H° d t s ° 2 is 
given by 



-Ef 



;PPe 



- 2J Z S Z - J+S- - S+J- 



AtA e + AtA , (2.10) 



As commented above, H° d t so1 does not appear for the 
zigzag cases due to a momentum mis-match. Further- 
more, for a thin CNT with zigzag edges, H° d t s ° 3 is also 
vanishing and hence the effective Hamiltonian is simpli- 
fied. 

The bulk part of a metallic CNT is described 
by a simple tight-binding Hamiltonian H(±ko) = 

—t Yli^i [ c ± i c ± i+i + h.c] after the gauge transforma- 
tion. In describing the gapless excitations, we replace 
lattice fermion operators c-t^by slowly varying contin- 
uum operators ipL/R,±{x) as 



c±,i 



D +ik F x n 



4>l,±{x)+z FX ^R,±{x) , (2-13) 



where a x := v / 3ao/2, and x := 2X . Due to the open 
boundary conditions, however, the left and right movers 
are not independent. They satisfy a constraint 



1pR,±( x ) = -i>L,±(- X )> 



(2.14) 



which allows us to concentrate on only the left-moving 
sector, say. Focusing on the two gapless modes, the ki- 
netic part of the system is written by the slowly varying 
variables as 



r + L 



(2.15) 



where L — N x a x and vp = \fitao/2 is the Fermi veloc- 
ity. Here, the original system defined for x € [0, +L] is 
extended to x € [-L, +L] by the constraint, Eq. (|2.14|) . 
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Upon taking continuum limit, we are also allowed to 
set X = x = for the interactions that include bound- 
ary states with some renormalizations of couplings, since 
edge modes are exponentially localized at the boundary. 
The lattice fermion operators at the boundary are re- 
placed with the left moving continuum field as 



c±,i=i 



ipL,±(a x ) + e 



/a x ~2i s\xi(k F a x ) ?Al,±(0). 



(2.16) 



Correspondingly, operators made of lattice fermions 
J,p,AJ,A appearing in ff° n d t E ° 2 are replaced with con- 
tinuum counterparts as 



3 = ^L s -^-vwi 

a,ss' 



(2.17) 



with suitable renormalizations for couplings. At this 
level, also, differences between the Hubbard and the 
Coulomb interaction are irrelevant, since we keep the 
same terms for the both types of interactions. 



Effects of the bulk interactions: Open boundary 
bosonization 



Although we are interested in physics at boundaries, 
effects of electron correlation in the bulk regime -ff^"" 1 
should also be taken into account, since the bulk interac- 
tion is known to affect the scaling dimension of operators 
inserted at a boundary in ID correlated systems. The 
bulk interaction is well incorporated by the bosonization 
technique for ID systems. Following Ref. I28L we bosonize 
the theory with open boundary condition. Open bound- 
ary bosonization was previously used to discuss the cor- 
relation effects of armchair CNTs with boundaries, where 
there are no boundary states though pIliSLiffii 

Upon bosonization, electron operators are expressed in 
terms of scalar bosonic operators as 



4> as (x) = 



i 



--Vo 



(x), (2.18) 



where rj as is a Klein factor, and : • • • : denotes normal 
ordering. [Since it is enough to focus on the holomorphic 
sector, we simply write ^l„ s (x) = ^ as (x), (p Las (x) = 
tp as {x), etc., henceforth.] For convenience, it is better to 
introduce a new basis {f P ,±,fa,±} rather than ip± ^/i, 

Vasix) := ^-{(p p+ (x) + sip a+ (x) 

+aip p -(x) +ascp a -{x)y (2.19) 

When the effects of forward scatterings are taken into 
account, left and right movers are mixed through the 
Bogoliubov transformation. Since there exists the con- 
straint between the left and right movers, bosonization 



rules take non-local form, where ip as in Eq. (|2.18l) is, 
after the Bogoliubov transformation, given by 

<Pjs(x) -> -\yK~ s [Vjs{+x) -<p jS {-x)} 

+ ^== [<Pji(+x) + Vjs{-x)] }, (2.20) 

where j = pja and 5 = +/—■ Especially at the boundary 

x = 0, 

(o) = ^ P+ (o)/V 

K p+ + scp a +(0)/ \/K a+ 

+a^ p _ (0)/ s/KpZ + asip a - (0)/ y/W^} . 

(2.21) 

Parameters Kj$ are the Luttinger parameters for each 
mode - that is, a coefficient of the Bogoliubov transfor- 
mation. The long-range bulk forward scattering strongly 
renormalizes the charge symmetric mode (p, +), and its 
Luttinger parameter is estimated to be K p+ ^0.2 for a 
CNT with the Coulomb interaction^ On the other hand, 
for the other modes, Kj$ is almost equal to unity if wc 
neglect the back scatterings and the umklapp scattering. 

The bulk interactions affect physics at boundaries 
through the modifications of the Luttinger parameters. 
More precisely, they alter the scaling dimensions of the 
operators inserted at boundaries. The scaling dimensions 
of J ± and At, A in Eq. (|2.17() . which we call x± and xa, 
respectively, are equal to unity in the absence of the in- 
teractions. However, they are now given by 



K 



a + 



k: 



K 



XA 



P+ 



k: 



(2.22) 



which are dependent on the Luttinger parameters. Wc 
see that superconducting pairing operators A^ , A are 
made to be strongly irrelevant by the strong bulk re- 
pulsive interactions. 

The bosonized expression for the bulk part of the 
Hamiltonian is 



TTbulk i rrbulk/ 

+ L 



rrbulk 



dx : Jjs(x) 2 :,(2.23) 



where Jjs = — dxfjs/^/^ represents a [/(l)-current and 
Vjs is a velocity for each collective mode. The velocity 
for the charge symmetric mode is strongly enhanced as 
v p+ ~ vf/K p+ , while the velocity is almost equal to the 
Fermi velocity for other modes, Vj$ ~ Vf, (J, S) ^ (p, +). 
H^ k ' represents the part that cannot be written as a 
current-current interaction, in the presence of which the 
modes cr± and p— are made to be gapped away from half 
filling, whereas all kinds of excitations are gapped at half 
filling due to the umklapp scattcrings>i£ii§*i£ 
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III. CASE OF ONE BOUNDARY STATE 



r c e , wc obtain a set of RG equations 



Having established low-energy effective theories that 
describe correlation effects at boundaries, we now turn 
to specific examples. We first consider a thin metallic 
CNT with only one boundary state, such as (6, —6) CNT 
with zigzag edges or (3, —3) CNT with bearded edges. 
Especially, we focus on the former example near half fill- 
ing as it is the simplest case in that its boundary state 
f-n(i) appearing for k y a$ = ir is completely localized at 
boundaries, f WlP (i) = £j,i£ p ,. [Eq. l|2.3[l ] . 

The most general expression for low-energy effective 
theories is now reduced to the following Hamiltonian: 



H 
Hi 



rrbulk I rrbulk/ 
H + H int 



bulk/ 
nt 3 



H°«r =: H + H! + H! 
dx : Jjs(x) 2 : +e e p e + U e n^ni, 

- L 



v F X p 



p e p + V F X Z S Z J Z 



s + j- + J+S~ 



(3.1) 



where J, p are defined in Eq. (|2.17|) . Superconduct- 
ing pairing operators A ' , A are dropped as they are 
strongly irrelevant. Also, S and p e arc reduced to 
& = 2~2 SS > e t "'ss'/ 2 e s , and p e = J2 s e t e s' respectively. 
An ultraviolet cutoff t c is introduced, which is esti- 
mated to be the inverse of the bandwidth, r c ~ 1/t 
{vft c ~ a x ). A p , Zi j_ represents a dimensionless cou- 
pling between the boundary states and conduction elec- 
trons. For a (6, —6) CNT with zigzag edges, ini- 
tial conditions for RG analysis (bare values of the 
couplings) are given by U e = U /N, e e ~ 0, and 
-v F X z /2 = -v F X±/2 = v F X p /2 = Aa x sin 2 (kFa x )U/N 
for the Hubbard interaction near half filling, whereas 
they are given by U e = V(0), e e ~ 0, —v F X z /2 = 
—v F X±/2 = 4:a x s'm 2 (k F a x )V(ko + it), and v F X p /2 = 
4a x sin (k F a x ) [2V(0) - V(k + tt)] for the Coulomb 
interaction, where V(q) := Ny 1 J2y j e ~ 19Vj V(yj) an d 
V(yj — yji) = V p=p ' =m -i-i' = i-jji. Thus, if we switch off 
the interaction between boundary states and conduction 
electrons (Hi), the ground state at the boundary at (or 
slightly above) half filling is the state where one of the 
two boundary states is occupied. 

To see what happens when we switch on couplings be- 
tween boundary and conduction electrons, we perform 
a perturbative RG analysis up to one-loop order— via 
a Coulomb gas representation of the partition function, 
with Hq being an unperturbed part of the Hamiltonian. 
We neglect H^ lkl for the time being, and consider the 
case where the system is described by TL liquid, i.e., 
temperature above the gaps induced by H^ t ml . This is 
a good approximation for the systems off the half filled 
condition. Infinitcsimally rescaling the ultraviolet cutoff, 



dX p 
dl 

dX ± 
~df 

dh e 
~~dT ' 

dhu 
~df 




{l-x±)X± + 



I 4 4 2 



(3.2) 



where dimensionless couplings arc defined as h e 
hu ■= T c U e , and we performed the rescaling 



2tt 



v P + 



Xp, X z 



Up 



2- 



A i -» — 



1 


v F 


2if .+ 


V F 




2~H 






. v <y- _ 





A 



(3.3) 



We can treat A, and 



p non-pcrturbatively in the manner of Schotte and 



to simplify the RG equations 
A_ 

SchottoS, which actually gives an identical result to the 
above one-loop calculation for dh e /dl and dhjj/dl. How- 
ever, we prefer to treat X z , X p and A^ on an equal footing 
here. 

From the RG equations, we see that interactions in 
the spin sector A Zj ^ arc renormalizcd to zero as it is well 
known that the ferromagnetic, isotropic Kondo interac- 
tions are vanishing in the infrared. The bulk electron 
correlations are almost irrelevant for the RG equations 
of the spin sector. On the other hand, the repulsive bulk 
interactions profoundly affect the charge sector; they sup- 
press the charge fluctuations at boundaries. First, as al- 
ready commented, the scaling dimensions of A^ and A, 
which allow a pairwisc hopping from the bulk part to 
boundary states and vice versa, are made strongly ir- 
relevant. Second, they modify the values of couplings 
through the Bogoliubov transformation as seen from Eq. 
I|3.3(l . When the bulk repulsive interaction is very strong 
K p+ « 1, the bare value of A p is drastically reduced 
after the rescaling (|3.3f) . which amounts to U e — > +oo 
and U e + 2e e — > +oo as I — > +oo. Then, doubly oc- 
cupying a boundary state is prohibited in the infrared 
limit. Note also that even though we start from slightly 
below half filling h e > 0, for which no boundary states 
are occupied if we switch off the coupling between bulk 
and boundary electrons, couplings X p . z .± rcnormalize h e 
to — oo and hence one of two boundary states is occu- 
pied in the ground state, with the total spin carried by 
the boundary states being equal to 1/2. We numerically 
solve the RG equations (|3.2|) for the Coulomb interac- 
tion and confirmed that this is the case for the Coulomb 
interaction^. 

At the fixed point X z .± — > 0, e e — > — oo, and U e , U e + 
2e e — > +oo, the conduction electron is described as iso- 
lated TL liquid with the open boundary condition, where 
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tunneling density of states or Freidcl oscillation can ex- 
hibit characteristic behavio»i§i22*22i3i. Correlation effects 
between bulk and boundary states around this fixed point 
can be taken into account in a perturbative way. 

At very low temperature, H^. M causes energy gaps in 
the modes <r± and p— away from half filling, whereas all 
kinds of excitations arc gapped for half filling. Effects of 
iJ^" lk/ on boundary physics are twofold. First, although 
we have treated the Luttinger parameters as a fixed con- 
stant, they undergo renormalizations when the effects of 
H^_ ml are taken into account. The Luttinger parameters 
K a+ , K p -, and A'^j 1 always renormalize to zero, while 
K p+ renormalizes to zero only for half filling>i£ii§ii£ This 
effect is easily accounted for in the RG equations Ij3.2|l . 
which does not affect the conclusions as A* and A stay 
irrelevant. Second, H ] ™ lkl might generate extra operators 
at the boundary upon RG transformations. Then, RG 
equations should be traced with these extra operators. 
This is, however, beyond the present discussion. 

To get further insights, it might be better to adopt 
a complementary starting point rather than the above 
TL liquid model for higher temperature. In the lower- 
temperature limit, slightly away from half filling, a su- 
perconducting ground state on the honeycomb lattice can 
be formed, which originates from a rung-singlet state in 
the effective two-leg ladder model. Emergence of iso- 
lated states at boundaries can then be determined based 
on this specific pattern of singlet pairs, as in the valence- 
bond solid states in spin systems. 



IV. CASE OF TWO BOUNDARY STATES 



boundary near half filling is found to be the state with 
the total spin equal to unity. 

To see the effects of Hj, RG equations for the cou- 
plings are obtained in the same way as the case of one 
boundary state. Again, we focus on TL liquid regime 
and neglect H^ t lk ' . RG equations for A z , Aj_, h e , and hu 
are identical to those in the case of one boundary state. 
Then, A z and Aj_ become vanishing in the infrared limit, 
since they are initially ferromagnetic and isotropic. In 
addition, charge fluctuations are suppressed when there 
are the strong repulsive interactions in the bulk, since 
e e — > — oo and U e ,U e + 2e e — > +00, and hence doubly 
occupying a boundary is unfavorable. RG equations for 
K z .± and /, which determine the total spin carried by 
the ground state of e fermions, are given by 



dhl -h 9X 2 



dh 



(11 



fix* - A, 



dh 



dl 



Ai, (4.2) 



where Hk, ± '■= t c K z ^± and hj := t c I. We see that 
the Kondo couplings \ z ,± renormalize the exchange in- 
teractions between boundary states AT Zi j_, making it fer- 
romagnetic. Then, the ground state of the boundary 
states is polarized with the total spin equal to unity. 
When we include At, A in the RG analysis, they also 
give rise to a contribution in the perturbative expansion. 
In contrast to A Zj _l, they suppress the ferromagnetic cou- 
pling between the boundary states K Zi ±. However, the 
bulk interaction makes A T , A strongly irrelevant as stated 
above, and hence they do not affect the RG flow. 



As a next step, we consider a thicker metallic CNT 
with two boundary states: (9, —9) CNT with zigzag 
edges, for which boundary states appear for k y ao = 
— 87r/9 and +8tt/9. (Fig. Our main interest here is on 
whether or not the total spin carried by the two bound- 
ary states is non-zero. The effective Hamiltonian for this 
case is given by 



Hn = 



dx : Jjg(x) 2 : 



it 

+ Lp lP2 + K z S z 1 S% + ^ 



v F X p 



4 
i>fA 



nii nil + n 2|^2i 
p e p + v F \ z S z J : 



-{v F r c f 



s + j- 



S x S 2 H 
e e p e , 

j+s- 



(4.1) 



where we simply write ei := e +s „/g, e2 := e_ 87r / 9 . Again, 
superconducting pairing operators A T ,A are dropped as 
they are irrelevant. Initial conditions are given by K z = 

K x = -I = -2U e < 0, e e - 0, and A z = A_l A p < 0. 

Then, if we neglect the couplings between conduction 
electrons and boundary states, the ground state for the 



V. CONCLUSION 

To conclude, we have investigated correlation effects of 
(N, —N) CNT's with boundaries. Low-energy effective 
Hamiltonians were established by taking into account 
boundary states at the fermi energy. Due to the special 
nature of the boundary states, the differences between 
the Coulomb and the Hubbard interaction are found to 
be small. We then discussed specific examples where only 
one or two boundary states appear. In the infrared limit 
in RG analyses, Kondo-like couplings between conduc- 
tion electrons and boundary states are shown to be van- 
ishing, and the bulk conduction electrons and boundary 
states are completely decoupled. Doubly occupying a 
boundary state becomes unfavorable in the infrared near 
half filling since the strong repulsive bulk interactions 
renormalizes the interactions at the edge. As the bound- 
ary states do not dissipate through the coupling with 
conduction electrons, they can be directly observed by 
local probes such as scanning tunneling microscope. The 
boundary states also manifest themselves in transport 
experiments as the conduction electrons and the bound- 
ary states give rise to independent contributions to, say, 
tunneling density of states. 
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Furthermore, the ground state at the boundary is a 
highest spin state S = 1/2 for the case of one boundary 
state and S = 1 for two boundary states, since ferro- 
magnetic couplings between the boundary states are en- 
hanced by the interactions between conduction electrons 
and the boundary states. Then, the boundary states, as 
a whole, behave as a localized moment which gives rise 
to a Curie- Weiss-like contribution to the magnetic sus- 
ceptibility, apart from the bulk conduction electrons. 

The result obtained here is consistent with spin polar- 
ization found in DMRG study for a thin semi-conducting 
CNT—, and mean-field theory^ or a LSDA calculation^ 
for 2D sheet geometry. 

At lower temperature, the formation of a spin-gapped 
ground state is suggested by the previous studies. The 
fate of the spin polarization, found here for the temper- 
ature above the gaps, is left as an open question. 

Finally, we comment on an extra gapless mode other 
than the modes for k y ao = ±2ir/3 treated in the present 



paper. When one considers the hybridization between 
(7 — 7r orbitals, which the simple tight-binding approx- 
imation adopted here does not correctly capture, there 
can appear a gapless mode for a very thin CNT, as sug- 
gested by band calculations^— In a realistic situation, 
effects of the extra gapless mode should be taken into ac- 
count, which can be possible along the lines of the present 
discussions. 
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